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Abstract 

LISA (Laser Interferometer Space Antenna) is a mission to detect and study low- frequency cosmic 
gravitational radiation through its influence on the phases of six modulated laser beams exchanged 
between three remote spacecraft. We previously showed how the measurements of some eighteen 
time series of relative frequency or phase shifts could be combined (1) to cancel the phase noise 
of the lasers, (2) to cancel the Doppler fluctuations due to non-inertial motions of the six optical 
benches, and (3) to remove the phase noise of the onboard reference oscillators required to track the 
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photodetector fringes, all the while preserving signals from passing gravitational waves. Here we 
analyze the effect of the additional noise due to the optical modulators used for removing the phase 
(T) ' fluctuations of the onboard reference oscillators. We use the recently measured noise spectrum of 

an individual modulator (Klipstein et al. Ill) to quantify the contribution of modulator noise to 
. the first and second-generation Time-Delay Interferometric (TDI) combinations as a function of 

o. 

the modulation frequency. We show that modulator noise can be made smaller than the expected 



proof-mass acceleration and optical-path noises if the modulation frequencies are larger than ~ 682 



^ ■ MHz in the case of the unequal-arm Michelson TDI combination X\, ~ 1.08 GHz for the Sagnac 

b : 

TDI combination ct\, and ~ 706 MHz for the symmetrical Sagnac TDI combination Q\. These 
modulation frequencies are substantially smaller than previously estimated and may lead to less 
stringent requirements on the LISA's oscillator noise calibration subsystem. The measurements 
in [l| were performed in a laboratory experiment for a range of modulation frequencies, but we 
emphasize that, for the reference oscillator noise calibration algorithm to work, the modulation 
frequencies must be equal to the frequencies of the reference oscillators. 
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I. INTRODUCTION 



LISA (Laser Interferometer Space Antenna) is a three-spacecraft deep space mission, 
jointly proposed to the National Aeronautics and Space Administration (NASA) and the 
European Space Agency (ESA). It will detect and study low-frequency cosmic gravitational 
radiation by observing frequency shifts of laser beams interchanged between drag-free space- 
craft y. 

Modeling each spacecraft with two optical benches, carrying independent lasers, beam 
splitters and photodetectors, we previously analyzed the measured eighteen time series of 
frequency shifts (six one-way laser carrier beams between spacecraft pairs, six between the 
two optical benches on each of the three spacecraft, and six more by over- imposing phase 
modulations on the laser carrier beams between spacecraft pairs to monitor rates of the 
reference oscillators on each spacecraft). We showed that there exist several combinations of 
these eighteen observable which cancel the otherwise overwhelming phase noise of the lasers, 
and the phase fluctuations due to the non-inertial motions of the six optical benches, and 
also allow the removal of the phase noise from the onboard reference oscillators (or USOs - 
"Ultra-Stable Oscillators") required to track the photodetector fringes, while leaving signals 
due to passing gravitational waves 

The analysis in our previous work assumed that noise due to the electro-optic modulators 
used to implement the onboard reference oscillator noise calibration algorithm was negligible. 
Here we amend those results to a more realistic LISA operational configuration where the 
effects of the phase fluctuations of the modulators are explicitly included in the various time- 
delay interferometric (TDI, [4() combinations. An experimental investigation for estimating 
the magnitude of the phase noise introduced by a commercially available electro-optical 
modulator (EOM) on the phase measurement was recently performed by Klipstein et al. [l| ■ 
They showed that a modulation frequency of about 8 GHz would be adequate to suppress 
EOM noise in a one-way phase measurement to smaller than a budgeted noise level. Their 
conclusion was conservative in that it did not take account of the transfer functions of 
modulator noise to the TDI combinations. Here we derive the EOM noise transfer functions 
and show that smaller modulation frequencies than those previously identified could be 
used, implying less stringent phase noise requirements on the EOMs and perhaps simpler 
subsystem design. 
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In Section [IT] we give a brief summary of TDI and its implementation for LISA, and 
show that the expressions derived in for removing the Ultra-Stable Oscillator noise from 
the TDI combinations, valid for a stationary LISA configuration (so called "first-generation 
TDI"), can easily be generalized to the "second-generation" TDI combinations (i.e. those 
accounting for the the motions of the three spacecraft with respect to each other and around 
the Sun). 

In Section II III we derive the transfer functions of the modulator noises to the second- 
generation TDI combinations X\ (unequal-arm Michelson), ot\ (Sagnac), and d (sym- 
metrized Sagnac). Using the measurements in ij] on a particular EOM, we show that 
selecting modulation frequencies greater than m 682 MHz (for Xi), ~ 1.08 GHz (for «i), 
and ~ 706 MHz (for d) results in the power spectral density of the modulator noise being 
smaller than the optical-path and proof-mass noises of these TDI combinations. 



II. TIME-DELAY INTERFEROMETRY 



Equal-arm Michelson interferometer detectors of gravitational waves can observe gravita- 
tional radiation by canceling the much larger frequency fluctuations of the laser light injected 
into their arms. This is done by comparing phases of split beams propagated along the equal 
(but non-parallel) arms of the detector. The laser frequency fluctuations affecting the two 
beams experience the same delay within the two equal-length arms and cancel out at the 
photodetector where relative phases are measured. In this way gravitational wave signals 
of dimensionless amplitude less than 10~ 22 can be observed using lasers whose frequency 
stability can be as large as 10 -13 . 

If the arms of the interferometer have different lengths, as will inevitably be the case for 
space-based detectors such as LISA, simple differencing of the phases on the photodetector 
does not exactly cancel the laser phase fluctuations, p(t). The larger the difference between 
the two arms, the larger will be the magnitude of the laser phase fluctuations affecting the 
detector response. If L\ and L 2 are the lengths of the two arms, the laser relative phase 
fluctuations remaining in the response is equal to (units in which the speed of light c = 1) 

Ap(t) = p(t-2L 1 )-p(t-2L 2 ) . (1) 

In the case of LISA, whose lasers are expected to display relative frequency fluctuations 
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equal to about 10~ 13 / V H z in the milliHertz band, and whose arms will differ by a few 
percent [2j], equation ([1]) implies the following expression for the amplitude of the Fourier 
components of the uncanceled laser frequency fluctuations (an over imposed tilde denotes 
the operation of Fourier transform) 

|A^M|~|pH| 2w|(L 1 -L 2 )| . (2) 

At u/2n = 1CT 3 Hz, for instance, and assuming \Li — L 2 \ — 0.5 sec, the uncanceled fluctu- 
ations from the laser are equal to 6.3 x 1CT 16 / vHz. Since the LISA sensitivity goal is about 
10~ 20 /v / Hz in this part of the frequency band, it is clear that an alternative experimental 
approach for canceling the laser frequency fluctuations is needed. 

The solution of this problem is achieved by first "de-coupling" the two arms with the 
implementation of a multi-photo-receiver design in which, at each optical bench, the phase 
difference between the light entering the arm and the one exiting it is measured, time-tagged, 
and digitally recorded. By then properly time-shifting and linearly combining these phase 
measurements exact cancellation of the laser phase noise is again achieved. The sensitivity 
to gravitational radiation for the practical LISA case is essentially equivalent to that of an 
equal-arm Michelson detector. For a physical and historical description of this technique, 



called Time-Delay Interferometry (TDI), the reader is referred to [4| and references therein. 

In the case of LISA there are six beams exchanged between the three spacecraft, with the 
six relative phase measurements s»j (i, j = 1,2,3) recorded when each received beam is mixed 
with light from an independent laser on the receiving optical bench (that laser also being 
used to transmit light back along the same arm to the distant spacecraft 1 ). The frequency 
fluctuations from all the six lasers, which enter in each of the six Doppler measurements, 
must be removed to levels smaller than those of the secondary (proof mass and optical path) 
noises jjj in order to detect and study gravitational radiation at the predicted amplitudes 
(see Figured] for a description of the LISA geometry). 

As the three spacecraft forming the LISA array orbit the Sun, their systematic relative 
motions result in Doppler shifts of the laser frequencies. These offsets may be as large as 
about ±20 MHz and are measured at the receiving spacecraft via heterodyning. To limit the 



1 This is the simplest configuration; for technical reasons, schemes in which two or more of the lasers are 
remotely phase-locked are being considered by the LISA project. TDI works as well for these alternate 
measurement approaches also [5[. 
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FIG. 1: Schematic LISA configuration. Each spacecraft is equidistant from the point O, in the plane 
of the spacecraft. Unit vectors hi point between spacecraft pairs with the indicated orientation. 
At each vertex spacecraft there are two optical benches (denoted 1, 1*, etc.), as indicated. 

offset frequencies to this range either the lasers must be referenced to an atomic line (such 
as molecular iodine stabilization jfl) or a sophisticated locking scheme must be employed 

Q. 

A heterodyne measurement is "base-banded" by using a properly selected tracking fre- 
quency generated by a reference oscillator (the USO) using servo loop feedback. Although 
this procedure allows us to accurately track the phases of the photodetector fringes, it intro- 
duces into the phase measurements a residual noise due to the USO itself Q|. This USO noise 
is not negligible and a technique for removing it from the TDI combinations was devised (see 
j^l for more details). This technique requires the modulation of the laser beams exchanged 
by the spacecraft, and the further measurement of six more inter-spacecraft relative phases 
w comparing the sidebands of the received beam against sidebands of the transmitted beam 
I 

The time-keeping of the three spacecraft must be synchronized in an inertial frame, usu- 
ally taken to be the solar system barycentric frame. On each link, the up- and downlink delay 
times used in the TDI combinations thus differ due to large relativistic aberrations. This is 
first accommodated in the analysis by so called "modified generation" TDI combinations, 
and in "second generation" TDI where more general "flexing" of the LISA configuration is 
allowed. Following the arms are labeled with single numbers given by the opposite space- 
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FIG. 2: Schematic diagram of the LISA configuration involving six laser beams. Optical path 
delays taken in the counter-clockwise sense are denoted with a prime, while unprimed delays are 
in the clockwise sense. 

craft; e.g., arm 2 (or 2') is opposite spacecraft 2, where primed delays are used to distinguish 
light-times taken in the counter-clockwise sense and unprimed delays for the clockwise light 
times (see Figure (j2J)). Also the following labeling convention of the relative phase data will 
be used. Explicitly: S23 is the one-way phase shift measured at spacecraft 3, coming from 
spacecraft 2, along arm 1. Similarly, S32 is the phase shift measured on arrival at spacecraft 
2 alone arm 1' of a signal transmitted from spacecraft 3. Due to aberration and the relative 
motion between spacecraft, L\ 7^ L\ in general. As in |3j, we denote six further data streams 
Tij = 1, 2, 3), as the intra-spacecraft metrology data used to monitor the motion of the 
two optical benches and the relative phase fluctuations of the two lasers on each of the three 
spacecraft. 

The frequency fluctuations introduced by the lasers, by the optical benches, by the proof 
masses, by the fiber optics, by the USOs, and by the measurements themselves at the 
photo-detectors (i.e. the shot-noise fluctuations) enter the Doppler observables s^, with 
specific time signatures; see (3], Q] for a detailed discussion. The contribution s^ w due to 
GW signals was derived in in the case of a stationary array, and further extended to the 
realistic configuration [10] of the LISA array orbiting the Sun. 

The four data streams recorded at spacecraft 1, including Doppler effects, independent 
lasers, gravitational wave signals, optical path noises, proof-mass and bench noises, fiber 
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optics and USO phase fluctuations, are now given by the following expressions 
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For all the down conversions at spacecraft 1, the USO-generated frequency f\ is used, and 
it is multiplied by coefficients a 2 i, a 3i , c 2 i, and c 3 i that are the result of phase- lock loops 
driving numerically controlled oscillators to track the large (~ 20 MHz center frequency) 
offsets in the heterodyne phase measurements (the "beat notes") ^]. Thus the values of 
these coefficients are determined by the following expressions 

v 3 (1 - L 2 ) - v\ 
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(7) 



'1-L 3 ,) 



a 2 i 



fi 



(8) 



c 2 i 



-C31 



fl 



(9) 



Eight other relations, for the readouts at vertices 2 and 3, are given by cyclic permutation 
of the indices in equations ©-(JH))- 

Recent experimental results have indicated that onboard lasers can be effectively stabi- 
lized by referencing their frequency to that of molecular-line transitions such as those defined 
by molecular iodine (I 2 ). Besides providing a high reference frequency stability for the LISA 
onboard lasers, I 2 stabilization has the advantage of making the center frequencies of the 
lasers onboard each spacecraft essentially equal [6]. In the following we will assume the 
onboard lasers to be iodine-stabilized, and we will correspondingly set the coefficients c^- 
defined in equation to be identically equal to zero. 

The gravitational wave phase signal components in equations ([3]) and ([5]) are given by in- 
tegrating with respect to time the equations (1), (2) of reference |9J, which related frequency 
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shifts to metric perturbations. It is these that LISA is designed to measure at levels set by 
the optical path phase noise contributions, s° pt ' path , due mainly to shot noise from the low 
signal-to-noise ratio (SNR) in the links connecting the distant spacecraft. The measure- 
ments will be made with high SNR so that for them the shot noise is negligible. The other 
unavoidable secondary phase noise is due to proof mass non-gravitational perturbations, 
described here by the random displacement vectors A« and A*. 

The expressions of the phase measurements given by equations fl3]-[6]), when substituted 
into the second generation laser-noise-free combinations yield data that, although free of 
laser and motional phase noises, are now affected additionally by the USO phase fluctuations, 
which have been denoted {i = 1,2,3) in equations (J31 - ED • For instance, with a state-of- 
the-art USO displaying a frequency stability of about 1.0 x 10~ 13 in the milliHertz frequency 
band, the corresponding relative frequency fluctuations, (ji/i/i, introduced in the laser-noise- 
free data combinations would be equal to about 3.0 x 10~ 20 , several orders of magnitude 
above the secondary noises and LISA sensitivity goals [2|. 

The expressions of the gravitational wave signal, the USO and the secondary noise sources 
entering into X\ will in general be different from those entering into X, the corresponding 
"first generation" unequal-arm Michelson observable derived under the assumption of a 
stationary LISA array js), Q]. However, the magnitude of the corrections introduced by the 
motion of the array on the gravitational wave signal, the USO and the secondary noise 
sources entering into X% will all be proportional to the product of their time derivatives, the 
spacecraft relative velocities, and the LISA arm-length. At 1 Hz, for instance, the larger 
corrections (due to aberration) will be about five orders of magnitude smaller than the main 
terms. Since the amplitude of these corrections scale linearly with the Fourier frequency, 
we can completely disregard the time- and direction-dependence of the time-delays entering 
into these noise sources over the entire LISA band 

These considerations imply that the second generation TDI expressions for the grav- 
itational wave signal, the USOs and the secondary noises can be expressed in terms of 
the corresponding first generation TDIs. For instance, the second generation unequal-arm 
Michelson combination, (where the g-index indicates the inclusion of the USO noises) 
can be written in terms of the corresponding first generation unequal-arm Michelson com- 
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bination, X q (t), in the following manner [U 

X lq (t) = X q (t) - X q (t - 2L 2 - 2L 3 ) , (10) 

in which the time-dependence of the light-travel times can be completely disregarded. Equa- 
tion (fTUj) implies that the USO calibration procedure for the second generation TDI com- 
binations can actually be evaluated by simply considering the corresponding calibration 
expressions derived for the first generation TDI expressions given in 3J. For this reason, 
from now on we will focus our attention on the first generation combinations. 



III. MAGNITUDE OF THE MODULATOR NOISE INTO THE TDI COMBINA- 
TIONS 

n 

In the USO calibration scheme first proposed by Bender et al. [2], a second frequency 
is superimposed on each of the six main laser beams: specifically, beams originating at 
spacecraft i are modulated at the frequency /, generated by its onboard USO. The main 
carrier signal, together with a side-band (of intensity perhaps ten times lower than the 
intensity in the carrier [2]) are transmitted, and at the receiving spacecraft j they are 
heterodyned at a photo detector with the local laser beam uj also carrying side-bands 3j. If 
the frequencies of the sidebands are carefully selected to be larger than the laser frequency 
offsets (but to differ from each other by an amount smaller than the operational bandwidth 
of the photo detector) then the lowest two difference frequencies (or difference phase trains) 
can be distinguished and measured at the photo detector. These two difference phase time 
series are given respectively by the difference between incoming and outgoing carriers, and 
by the difference of the phases of their side-bands respectively. They are then independently 
further tracked with coefficients ay and fry [l^|. This process provides six additional 
data records, s[y Consider, for instance, the phase difference between the second signal 
transmitted by bench 3 and the second at the receiving bench 1* 



s 31 



[U Z + h) (1 - L 2 ) - V\ - /l - 6 2 l /; 



t + P3,2 ~ P*i 



+ <?3,2 - (1 + &21) qi + m 3i2 - ml (11) 

where, for simplicity, we have omitted the terms associated with the optical bench noises, 
the optical path and proof-mass noises, and the contribution from a possibly present grav- 
itational wave signal. Note that now the expressions for the measurements s„- include also 
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the phase noises mj,m* 



1, 2, 3 introduced by the modulators on the measured phase 



3| and here we amend that analysis. 



differences s'^. These terms were neglected in 

Note that the numerical coefficient 621, determined by the following equation 

fo + h) (i-l 2 )-K + /i) 



h 



(12) 



is distinct from a 2 i given by equation (J7J) (although they will be close if all the fi are close). 

The quantities r 2 ™ = (s 2 i — s' 21 )/ f 3 and similarly = (s 3 \ — s' 31 )/ f 2 (and cyclic 
permutations of their indices) enter into the algorithm presented in {3] for removing the 
USO noises from the various TDI combinations. Now they explicitly show their dependence 
on the modulator noises rrii, m* i — 1, 2, 3, and are related to the expressions introduced 
in [3]] (which did not include the modulator noises) by the following expressions 



A") 
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with the others obtained as usual via permutation of the spacecraft indices. 



A. Modulator Noise in the X-combination 



As a result of the presence of modulator noises in the s'^ measurements (and consequently 
in the combinations ) the first-generation TDI unequal-arm Michelson combination given 
in equation (27) of [3[ (which now will include combinations of the measurements) will 
be affected by the modulator noises. If we denote with the resulting new unequal-arm 
Michelson combination, it is easy to see that this can be written as the sum of X (as given 

n 

in [3[]) with terms due to the modulator noises. The resulting expression is equal to 



X^ = X - 012 h 
~m 

+ a 13 / 3 
- a 2 i /1 



"^1,22 



m. 



2,322 



+ 



m 



2.3 



mi 



h 



m\ - m 3i2 ^"3,2 - m 122 



h 



1,33 - m 3,233 m 3;2 - ml mi_ 



m 



2.3 



h 
m 3 o 



h 



+ 



/a 

/l 



/2 

+ «31 fl 



+ 



2,3 



fi 



fi 



mi 



m 



2.3 



h 



+ 



2,3 — m l,33 
li . 



(14) 
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The expression above leads to the following estimation of the spectral density of the noise, 
S X (.™) , in the TDI combination 

S x(m ) = S x + 4sin 2 (cjL) [S mi (a 12 + a 13 + a 31 ) 2 + S m * (a 12 + a 13 + a 2 i) 2 

-\-S m * a 12 + S m3 a 13 ] , (15) 

which has been written as the sum of the secondary noise spectra Sx of the combination X 
(not containing the modulator noises) and the spectra of the modulator noises themselves. 
Note that for a figure of merit equation ffTBl) conveniently assumes all USOs to have equal 
frequency /, and all the armlengths to be equal to a common value L. 

Assuming a maximum beat-note frequency offset of 20 MHz in each arm, the coefficients 
ciij can be taken all equal to a = 20 MHz//, where / is the USO or modulation frequency 2 . 
Given the measured modulator noise spectrum (figure 6 in [1]), from equation (|T5|) it is then 
possible to determine the value of the USO frequency which would make the modulator noise 
contribution smaller than the secondary noises (Q, Section IV) affecting the X combination 
over the entire LISA band. This can be done my maximizing over the angular frequency uj 
the following function (see equation [15] above) 



sin(cuL) ^/80 §4"4 x 20 MHz > ( 16 ) 

where S m (uj) is the spectrum of the noise associated with one modulator. Over the frequency 
band relevant to LISA, the measured S m (u) for a single EOM can be approximated by the 
following analytic expression (see figure 6 in l|) 

S m (u) = 2.8 x w" 1 cycles 2 Hz" 1 when: 10~ 4 < u/2n < 8 x 10~ 3 Hz 

= 5.5 x 10~ 8 cycles 2 Hz" 1 when : 8 x 10" 3 < u/2n < 3 x 10~ 2 Hz 

= 3.8 x 10~ 10 uj- 3 cycles 2 Hz" 1 when: 3 x 10~ 2 < u/2n < 5 x 10" 1 Hz 

= 1.2 x 10~ n cycles 2 Hz" 1 when : 5 x 10" 1 < u/2n < 1 Hz . (17) 

Note that the expressions of the power spectral densities of the noises entering into the X, 
a, and ( combinations are given in Section IV of [s| as fractional frequency fluctuations, 



2 In the experimental setup used by 1] to characterize an individual EOM's noise it was not necessary that 
the modulation frequency be the same as the USO frequency. We emphasize that in LISA's actual USO 
calibration subsystem the modulation frequency must equal the USO frequency; see Figure 5 of 
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log 10 (frequency, Hz) 

FIG. 3: Numerical comparison of modulation noise in X with a modulation frequency / equal to 
682 MHz (dashed-line) , 8 GHz (dash-dot line), and the combined secondary noises (solid line). 

and therefore need to be converted to the same units as S m before being used for direct 
comparison. 

After maximizing the function given in equation ( fl6l) we conclude that the modulation 
frequency / should be equal to or larger than 682 MHz for the modulator noise level in the 
X-combination to be smaller than the remaining secondary noises. A numerical comparison 
of modulation noise in X at this level, modulation noise if / = 8 GHz is adopted, and the 
combined secondary noises is shown in Fig. [3j 



B. Modulator Noise in the a-combination 



If we denote with the resulting new Sagnac combination, it is easy to see that this 



can be written as the sum of a (as given in 



and terms due to the modulator noises. The 



13 



resulting expression is equal to 

,M = a + [a 32 h + ai 3 / 3 ] 
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(18) 



The expression above leads to the following estimation of the spectral density of the noise, 
S a (m) , in the TDI combination 



S, 



(m) 



S a + (a 3 2 + 013) 2 (S mi + S m *) + (a 2 3 + CI12) 2 (S m * + S, 



"13/ 



+ |(ai3 + a a2 e^)| 2 ,S m2 + |(a 12 + ai3 e 4 ^)| 2 5, 



(19) 



which has been written as the sum of the noise spectra of the modulator noise-free combina- 
tion a and the contribution to the overall spectrum coming from the modulator noise. All 
the armlengths have been taken again to be equal to a common value L. 

Again taking the maximum beat-note frequency offset to be 20 MHz in each arm, the 
coefficients ajj can be treated as all equal to a = 20 MHz//. Now the function that needs 
to be maximized to determine the value of the USO frequency which makes the modulator 
noise smaller than the LISA secondary noises over the entire LISA band is equal to 



SJuj) 



(20 + 4cos(a;L)) x 20 MHz . 



(20) 



By calculating the maximum value of the above function over the LISA operational band we 
conclude that a modulation frequency / equal to or larger than « 1.08 GHz will suppress 
the modulator noise in the a-combination to a level smaller than the secondary noises. A 
numerical comparison of the resulting modulation noise in a with the secondary noises is 
shown in Fig. HI 



C. Modulator Noise in the ^-combination 



We denote with C (m) t 
the sum of £ (as given in 



re new symmetrized Sagnac combination, and write it again as 
3)]) with added terms due to the modulator noises. The resulting 
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FIG. 4: Numerical comparison of modulation noise in a with a modulation frequency / equal to 
1.08 GHz (dashed- line), 8 GHz (dash-dot line), and the combined secondary noises (solid line). 



expression is equal to 
C M = C+y(a 2 i-a 3i ; 
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In order to identify the minimum modulation frequency / that suppresses the modulator 
noise below the secondary noises, it is useful to set = ciji, which is allowed if there is no 
flexing and the six lasers are stabilized to the same I 2 transition frequency. Again assuming 
a maximum beat-note frequency offset of 20 MHz in each arm, the coefficients can be 
taken as a = 20 MHz//. A careful analysis shows that the function which needs to be 
evaluated to identify the value of the USO frequency / required to make the modulator 
noise smaller than the LISA secondary noises is 



'16 [2 - cos(2tt/L)] S m (u) 
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15 
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FIG. 5: Numerical comparison of modulation noise in £ with a modulation frequency / equal to 
706 MHz (dashed-line), 8 GHz (dash-dot line), and the combined secondary noises (solid line). 



By calculating the maximum value of the function above over the LISA band we conclude 
that a modulation frequency / equal to or larger than 706 MHz will suppress the modulator 
noise in ( to a level smaller than the secondary noises. A numerical comparison of the 
resulting modulation noise in ( is shown in Fig. 
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